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Abstract 

We discuss the relation between ensemble and time averages for quasistationary 
states of low-dimensional symplectic maps that present remarkable analogies with 
similar states detected in many-body long-range-interacting Hamiltonian systems. 
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1 Introduction 



A possible approach to complex phenomena, when complete solutions are not 
available, is to reduce the problem into a simpler one that still retains the im- 
portant features of interest, and to use the simplified version as a model that 
gives some insight about the original issue. Along these lines low- dimensional 
maps are useful tools, widely used in the literature, both to study specific 
behaviors of nonlinear systems with few degrees of freedom, and to under- 
stand the macroscopic properties of much larger systems. The macroscopic 
behavior that we are addressing in this paper is the emergence (depending 
on initial data) of quasistationary states (QSSs) characterized by an anoma- 
lous almost-constant value of macroscopic variables (like the average kinetic 
energy), before a crossover to equilibrium. These QSSs where for example dis- 
played in the dynamics of long-range iV-body Hamiltonian Mean Field (HMF) 
models where it has been shown the emergence of long-standing phases with 
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non-Gaussian velocity distributions and with an anomalous value of the tem- 
perature [1], before the Boltzmann-Gibbs (BG) equilibrium is attained. The 
physics associated with these QSSs is extremely rich, for example various con- 
nections with nonextensive statistical mechanics [1], with aging [2] and glassy 
dynamics [3] have been pointed out. 

Using a properly defined 'dynamical temperature' in coupled symplectic maps, 
it is possible to exhibit QSSs similar to those observed in the HMF model [4,5]. 
In fact, the presence of well known mechanisms (see, e.g., [6,7]) that are related 
to the Kolmogorov-Arnold-Moser (KAM) theory provides a possible theoreti- 
cal background for the explanation of such anomalous phases. Of course, the 
new interest in these results stems from the possibility of better understanding 
the dynamical mechanisms that lead to violations of the classical BG theory, 
such as the one discussed in [1] . In this paper we restrict to low-dimensions, ad- 
dressing the connection between ensemble and time averages. The former were 
already studied in [4]; the latter are obtained here following the dynamical be- 
havior of a single trajectory, in a similar way as a thermometer time-integrates 
the interactions with the system to provide the output temperature. 



2 QSSs using ensemble averages 



Symplectic maps are convenient tools for the study of Hamiltonian systems, 
since they are the result of a Poincare section in the phase space of a Hamilto- 
nian system. As pointed out in [4], QSSs similar to those detected in the HMF 
model can be reproduced in low-dimensional symplectic maps using ensemble 
averages. In this section we report the main results. 

In connection with the KAM theory, one of the most studied symplectic maps 
is the standard map [8], which is 2-dimensional. Since important qualitative 
changes in the topology of phase space occur when the Hamiltonian has more 
than two degrees of freedom (e.g., Arnold diffusion), it is also important to 
address this more general context by considering for example a 4-dimensional 
symplectic map obtained by coupling two standard maps: 



Pl {t + 1) = Pl (t) + § sin[27T0i(i)] (mod 1), 

(1) 



p 2 (t + 1) = p 2 (t) + § sm[2n9 2 (t)] (mod 1), 



0i(t + l) =p 1 (t+l) + e 1 (t) + bp 2 (t+l) (modi), 
2 (t + l) = p 2 (t + 1) + 2 (t) + b Pl (t + 1) (modi), 

where a±, a 2 , b e R, t — 0, 1, b is the coupling constant, 0, may be regarded 
as an angular variable, and as an angular momentum (i = 1,2). In the 
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Fig. 1. Case (A) 'Ensemble' QSSs in the standard map. (a) T e (t) for typical values 
of ai. We start with 'water bag' initial conditions (M = 2500 points in < 0\ < 1, 
pi = 0.5 ± 5 10~ 4 ). Inset: Inverse crossover time t c (inflection point between the 
QSS and the BG regime) vs. l/(ai — a c ) 2 ' 7 ■ (b) Time evolution of the ensemble 
in (a) for a\ = 1.1 (first row) and PDF of its angular momentum (second row). 
t = 0: initial conditions; t = t\ = 500: the ensemble is mostly restricted by cantori; 
t = ti = 10 5 : the ensemble is confined inside KAM-tori. See [4] for further details. 

following we will take into account two cases: (A) for b = we consider only 
coordinates 1 , that is, the 2-dimensional standard map; (B) we fix b = 2 and 
ai = 0,2 = a in order to analyze the 4- dimensional map (see [4] for details) . The 
standard map is integrable when a x = 0, while chaoticity rapidly increases with 
|ai|; we also remind that for positive a\ smaller than a c = 0.971635406... the 
(macroscopic) chaotic sea is disconnected, because of the presence of invariant 
trajectories, called KAM-tori, that span the whole interval 9 e [0, 1]. 

With some similarity with the so-called 'water bag' initial data that produce 
QSSs in the HMF model [1], we consider at t = a statistical ensemble of M 
copies of the map with arbitrary 0, and Pi randomly distributed inside small 
regions of the interval [0,1]. Usually, in the connection between statistical 
mechanics and dynamics for systems with diagonal kinetic matrix and zero 
average momentum, the temperature is set proportional to the average square 
momentum per particle (see, e.g., [9]). As we address situations with nonzero 
'bulk' motion, the analogous concept, here called dynamical temperature, can 
be denned as the (specific) variance of the total angular momentum: 

1 d/2 

^)^^E((P^))-(P^)> 2 ), (A):d = 2, (B):d = 4, (2) 



where () means ensemble average. The qualification dynamical is used since 
this definition purely descends from dynamics and not from a thermal contact 
with a thermometer. The temperature associated with the uniform distribu- 
tion in the entire phase space can be called BG temperature and it is given, 

for both cases, by T BG = ^ £?=i f/o d Vi v\ - (lo dpi Pif] = 1/12 ~ 0.083. 
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Fig. 2. Case (B) 'Ensemble' QSSs in two coupled standard maps, (a) T e (t) for b = 2 
and typical values of a. We start with 'water bag' initial conditions (M = 1296 
points with < 61,62 < 1, and p\,p2 = 0.25 ± 5 10~ 3 ). Inset: Inverse crossover 
time t c vs. 1/a 5 ' 2 . (b) Time evolution of the fractal dimension of a single initial 
ensemble in the same setup of (a). See [4] for further details. 

For large values of |ai| (case (A)) or \a\ (case (B)) the system is typically 
mixing and it rapidly relaxes the dynamical temperature to Tbg- Reduction 
of I ai| or \a\ causes the formation of islands, barriers and partial barriers in 
phase space. The latter are responsible for the appearance of the QSSs, since 
they confine the ensemble inside a limited volume during a certain time. If the 
initial data are defined inside a singly connected volume that does not contain 
islands and if the the available portion of phase space for the final relaxation 
of the ensemble is such that its projection over the (pi,p2)-plane (pi-axis in 
case (A)) is uniform, the final temperature coincides with Tbg- We remark 
that this is possible only if the dimension of the map is larger than two [4]. 
Fig. 1 and Fig. 2 display this behavior for cases (A) and (B) respectively. It 
is also important to notice that a nontrivial fractal dimension characterizes 
the QSSs, while, once the cross over to Tbg is performed, the dimension of the 
ensemble coincides with that of the phase space (see Fig. 2 (b)). This indicates 
a strong violation, during the QSSs, of the equal-a-priori probability postulate 
on which BG statistical mechanics relies. 



3 QSSs using time averages 



Given a dynamical function / that takes different values f(xo, t) along a single 
trajectory started in the phase space point xq, we define its time average 
f(xo,t, t ) at time t — t during the observation time t = tf — U as 

7KM ) = gq^^ , (3) 
it — is 
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Fig. 3. Case (A) 'Time' QSSs in the standard map. (a) Tt(i) for typical values of 
a\. Evolution for a single orbit started in xq = (0.3,0.5), with observation interval 
t Q = 10 2 . Dashed lines indicate Tbg- (b) Averages over 10 5 realizations of the kind 
in (a), with initial data inside < 6\ < 1: p\ = 0.5 ± 5 10~ 4 . 

where t« and tf are respectively the initial and final observation time, and 
t = t i + t /2. This definition mimics the action of a detector testing the 
behavior of the system over a portion of a particular trajectory. Within this 
context, the definition of the dynamical temperature is reproduced as 



1 d / 2 _ 

T t (x ,t,t o ) = — £(P 2 -P')> (A):d = 2, (B):d = 4. 



(4) 



For uniform distributions of momenta over the whole phase space during the 
observation interval t Q we obtain once again T t = T BG . 

When the system is sufficiently chaotic, as expected, T t oscillates around Tbg- 
On the other hand, if |ai| (case (A)) or \a\ (case (B)) are diminished, T t 
exhibits, instead of the two plateaux structure, an intermittent behavior. This 
indicates the presence of sequences of partial barriers along the phase space 
trajectory, that alternate rapid and slow diffusion processes. Notice that T t 
can jump to larger values than Tbg (i- e -j the variance of the time-distribution 
of total momentum can become larger than that of a uniform distribution). 
Fig. 3(a) and 4(a) show this behavior for cases (A) and (B) respectively, for 
typical orbits. 

Another interesting quantity is the average value of T t (xo,t,t ) calculated for 
different orbits x : T t av (t,t ) =(T t ) Xo . Of course, T t av coincides with T e for 
t Q = 1, while for to >> 1 it represents the average result of different mea- 
surements and it is a typical observable calculated in the ergodicity analysis 
(see. e.g., [9]). Fig. 3(b) and 5(a) display that in this way we obtain a be- 
havior qualitatively equivalent to the ensemble averages of Fig. 1(a) and Fig. 
2(a), but typically with lower temperatures. Interestingly enough, 5(b) seems 
to indicate a tendency for Tf v to approach T e , as t Q increases. 
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Fig. 4. Case (B). (a) T t (t) for b = 2 and typical values of a. Evolution for a single 
orbit started in xq = (0.5, 0.25, 0.8, 0.25), with observation interval t Q = 10 2 . Dashed 
lines indicate Tbg- (b) Phase space analysis of the trajectory in (a) for a = 0.3. Dots 
represent the projection of the orbit on the plane (6±,pi). 

4 Conclusions 

We have addressed a simple connection between nonlinear dynamical systems 
and thermostatistics. Using well known results in symplectic maps (see, e.g., 
[6,7]), we pointed out the presence of QSSs with remarkable similarities with 
analogous ones detected in the HMF model [1] . Common features include the 
presence of an anomalous temperature-plateau in correspondence of particu- 
lar values of a control parameter (ai or a for the maps, the specific energy 
for the HMF). QSSs appear in both cases for particular classes of initial data 
and after a certain amount of time a cross over to equilibrium is observed. As 
the number N of coupled elements goes to infinity, the duration of the QSS 
diverges (see [1] for the HMF and [5] for the maps), opening the possibility 
for this anomalous effect to be physically relevant if the thermodynamic limit 
is taken before the infinite-time limit. All these similarities endow with re- 
newed interest the analysis of dynamical mechanisms of symplectic maps. For 
example, an important result is that during the QSSs of low- dimensional sym- 
plectic maps the phase space occupation exhibits evidences of fractalization 
that disappear once equilibrium is attained. 

The 'core' of equilibrium statistical mechanics consists in the elimination of the 
time variable as can be justified for example by the ergodic hypothesis. In this 
paper, working with low-dimensional symplectic maps, we have illustrated that 
time and ensemble averages, due to the complexity of phase space, originate 
different behaviors. A single time average of the dynamical temperature, a 
variable that was introduced in [4] in order to enable a natural comparison 
with many-body Hamiltonian systems, displays an intermittent behavior. In 
contrast, averaging many of these histories with a large enough observation 
time we observed a two-plateaux structure qualitatively similar to the one 
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Fig. 5. Case (B) 'Time' QSSs in two coupled standard maps. T® v (t) obtained aver- 
aging 10 5 realizations of the kind in Fig. 4, with initial data inside < #i,#2 < L 
Pi,P2 = 0.25 ±5 10~ 3 . (a) The observation interval is t Q = 10 2 . (b) The observation 
interval is t Q = 10 3 . 

obtained with ensemble averages. 
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